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ABSTRACT 

There have been a number of attempts to measure the expansion rate of the uni- 
verse at high redshift using Luminous Red Galaxies (LRGs) as " chronometers" . The 
method generally assumes that stars in LRGs are all formed at the same time. In this 
paper, we quantify the uncertainties on the measurement of H(z) which arise when 
one considers more realistic, extended star formation histories. In selecting galaxies 
from the Millennium Simulation for this study, we show that using rest-frame criteria 
significantly improves the homogeneity of the sample and that H (z) can be recovered 
to within 3% at z ~ 0.42 even when extended star formation histories are considered. 
We demonstrate explicitly that using Single Stellar Populations to age-date galaxies 
from the semi-analytical simulations provides insufficient accuracy for this experiment 
but accurate ages are obtainable if the complex star formation histories extracted from 
the simulation are used. We note, however, that problems with SSP-fitting might be 
overestimated since the semi-analytical models tend to over predict the late-time star- 
formation in LRGs. Finally, we optimize an observational program to carry out this 
experiment. 
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1 INTRODUCTION 

An important constraint on cosmological models is the evo- 
lution of the Hubble parameter, which is defined as: 



H(z) = - 



dz 



dt 



(1) 



Most tests of cosmology use measures of the luminosity dis- 
tance or the angular diameter distance, which include an 
integral over H(z), but measuring H(z) directly can po- 
tentially provide a more direct constraint on cosmological 
parameters (Jimenez & Loeb 2002). H(z) can be deter- 
mined at high redshifts by measuring the time-interval, At, 
corresponding to a redshift interval, Az, to obtain an ap- 
proximation of the derivative in equation (1). A number 
of authors have attempted to use this method (Ferreras, 
Mechiorre, & Silk 2001, Jimenez & Loeb 2002, Jimenez et 
al. 2003, Capozziello et al. 2004, Ferreras, Melchoirre, & 
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Tocchini-Valentini 2003, Simon et al. 2005, Dantas et al. 
2007, Verkhodanov, Parijskij, & Starobinsky 2007, Samushia 
et al. 2009) to track the evolution of H(z) as a function of 
redshift, and place constraints on cosmological parameters. 
Most recently, Stern et al. (2010) extended the studies of 
Jimenez et al. (2003) to a larger, more homogeneous sample 
to make the measurement of H(z) from LRGs. 

Jimenez & Loeb (2002) suggested using spectroscopic 
age-dating of Luminous Red Galaxies (LRGs) to measure 
dz/dt. LRGs are galaxies selected from the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) by apparent magnitude to 
trace the evolution of a homogeneous, volume-limited sam- 
ple of red galaxies (Eisenstein et al. 2001). Jimenez et al. 
(2003) used LRGs from the SDSS to measure H a = 69 ± 12 
km s _1 Mpc -1 (the Hubble parameter at z=0). Simon et al. 
(2005) used the same method to find H(z) as a function of 
redshift for < z < 2. In both cases, the authors assume 
that LRGs are drawn from the same parent population with 
the bulk of their stars formed at high redshift and fit single- 
burst equivalent ages to the galaxies. Although there is evi- 
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dence that LRGs make up a fairly homogeneous population 
of galaxies consisting mostly of old populations, Roseboom 
et al. (2006) have shown that this is not completely true. 
Furthermore, Barber et al. (2007) find that these galaxies 
have a range of formation ages and stellar histories. Jimenez 
and Loeb (2002) note that galaxies would not have formed 
at exactly the same time and explore what effect this would 
have on one's ability to discriminate between different mod- 
els for w(z). They do not, however, discuss the consequences 
of star formation being extended over a period in each galaxy 
(they assume all stars form at exactly the same time in a 
galaxy). This extended star formation history (SFH) needs 
to be considered when one defines and measures an age for 
a galaxy. In Jimenez et al (2003), the question is addressed 
in a little more detail: they simulate LRG spectra and add 
a contribution from a small fraction of stars formed at late 
times; they then age-date the galaxies and conclude that the 
age "edge" is not significantly different from the theoretical 
expectation. They do not attempt to quantify what uncer- 
tainties arise as a result of the extended star formation. 

Given the importance of correct estimates of uncertain- 
ties on the measured H(z), we explore, in this paper, the un- 
certainties that arise from LRGs having extended and varied 
SFHs. We use LRGs identified in the Millennium Simulation 
(MS, Springel et al. 2005) since the only way to reliably de- 
termine uncertainties is in a simulation where the SFH and 
input cosmology are known. The MS is currently one of the 
most detailed simulations of the evolution of galaxies and 
large-scale structure and has been tuned to reflect observa- 
tional data. Almeida et al. (2008) studied LRGs in simula- 
tions and showed that the luminosity function and the clus- 
tering of LRGs are reasonably well-modelled in their semi- 
analytical models of galaxy evolution, which also reproduce 
many other observables. The parameters of the de Lucia et 
al. (2006) semi-analytical model have also been successfully 
tuned to match the luminosity, colour, and morphology of lo- 
cal elliptical galaxies, the more massive candidates of which 
would be LRGs. 

On the other hand, galaxies are not modelled perfectly 
in the simulations: there are indications that model galax- 
ies assemble later than observations may suggest (Collins et 
al. 2009) but more reliable estimates of stellar mass are re- 
quired. It should be noted that these models overestimate 
the small amount of evolution seen in the properties of LRGs 
(Wake et al. 2006, 2008) and likely over predict the amount 
of growth occurring in these sources (Masjedi et al. 2008). 
Specifically for LRGs, Barber et al. (2007) find similar star 
formation rates, but younger ages than the de Lucia models. 
In addition, properties of local galaxies can be reproduced 
with different models for galaxy formation (e.g. Keres et 
al. 2009). However, the simulations provide viable models 
of extended SFHs, allowing us to quantify associated un- 
certainties and optimize analysis methods before applying 
them to observational data. 

In this paper, we optimize the selection of LRGs based 
on their properties in the MS and examine how well H(z) 
can be recovered from their mass-weighted age. To explore 
and quantify the uncertainties associated with extended 
SFHs, we then investigate how well H (z) can be determined, 
if the mass- weighted age of these LRGs were measurable, us- 
ing three different methods. In practice, mass-weighted ages 
are not directly measurable and accurate age-dating using 



LRG spectra is potentially a bigger challenge for the cos- 
mic chronometer approach. We carry out an initial study 
of age-dating, and its implications for H(z) measurements, 
using the stellar population models of Bruzual & Chariot 
(BC2003, 2003) to synthesize spectra. In contrast to fitting 
only single stellar population models (SSPs) (e.g. Simon et 
al. 2005) or dual-burst population models (Jimenez et al. 
(2003)) to determine the LRG spectral ages, we make a 
comparison of SSP results to those using more complicated 
SFHs. In a companion paper (Olivier et al. 2009), we delve 
much deeper into issues of age-dating galaxy spectra. 

The outline of the paper is as follows: In 321 we discuss 
the extraction of LRGs from the MS. In Sj3]we investigate 
how well H(z) is recovered in the MS using LRGs identi- 
fied in the simulation. In §3J we describe the modelling and 
age-dating of LRG spectra. In ij5] we discuss an observa- 
tional program to measure H{z). We summarize our results 
and conclusions in SJS] Throughout this work, we adopt the 
cosmology used in de Lucia et al. (2006) of Q. m = 0.25, 
n A = 0.75, and h = 0.73. 



2 IDENTIFYING LRGS IN THE MILLENNIUM 
SIMULATION 

The Millennium Simulation is a dark-matter-only N-body 
simulation run using the Gadget code of Springel et al. 
(2005). The simulation has I0 10 particles distributed in a 
500/i _1 Mpc 3 box. Dark matter halos are identified in snap- 
shots of the simulation at 63 different times and merger trees 
of halos are constructed. Semi-analytical modelling of galaxy 
evolution is carried out using the merger trees, providing a 
catalogue of galaxies in each snapshot, where each galaxy 
has an array of observational properties associated with it 
(Baugh et al. 2005, Bower et al. 2006, de Lucia et al. 2006). 
Properties of the halos and simulated galaxies are stored in 
a database which can be accessed from the web. 

Here, we focus on the galaxy catalogues constructed us- 
ing the semi-analytical models of de Lucia et al. (2006). We 
identify simulated galaxies with LRG properties in various 
snapshots (corresponding to a particular set of redshifts near 
z ~ 0.5) and then extract the star formation histories. The 
model does reproduce many properties of local galaxies and 
LRGs (Croton et al. 2006, de Lucia et al. 2006), but a num- 
ber of weakness are apparent as outlined in the introduction. 

In addition to the de Lucia models, we examined the 
models by Bower et al. (2006, hereafter Durham models), 
which also uses the Millennium Simulation. However, the 
two models are different in the implementation of their semi- 
analytical models. For our purposes, the most relevant differ- 
ence between the models is in their handling of gas cooling, 
star formation, and feedback as these are directly related 
to the star formation histories of the galaxies. We save the 
comparison of other semi-analytical models and how they 
would affect the calculation of H(z) for future work. 

2.1 Selecting LRGs using SDSS criteria 

Initially we extracted LRGs using the apparent magnitude 
constraints for SDSS data from Eisenstein et al. (2001). Our 
target redshift is close to the redshift where the SDSS 'Cut 
I' must be replaced with the 'Cut IF colour selection criteria 
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(At z ~ 0.4, the Balmer break moves into the r-band and 
a different colour selection is required to identify LRGs). In 
Fig. [U we plot the average SFH for three redshift snapshots 
at z — 0.32, 0.46, 0.51. The SFHs, especially once normalized 
to the SFH at z = 0.32, show significant differences. Impor- 
tantly, there is no monotonic behaviour for successive z. In 
Table 1, we summarize the average formation parameters for 
the population. There is a significant dispersion in formation 
redshift (zf), the redshift of peak star formation (z p ), and 
the redshift of last major star formation epoch (z e ), which 
is defined as the redshift where the star formation rate falls 
below 10 Mq yr~ x . For galaxies in the MS at z — 0.51, al- 
most 30% have had a burst of star formation greater than 5 
M^yr -1 in the previous gigayear and up to 19% are greater 
than 3cr away from the average star formation history. 

In addition, the change in the average mass-weighted 
age, Age mm , does not correspond to the change in the age 
of the Universe at these redshifts. Between z — 0.32 — 0.46, 
the age of the simulated Universe changes by 1.13 Gyr, 
whereas the change in Age mw corresponds to 1.27 Gyr. For 
2 = 0.46 — 0.51, the values are 0.26 and 0.18 Gyr, respec- 
tively. A key assumption of this technique is that the galax- 
ies are passively evolving. If true, then the change in their 
age with redshift should corresponding to the change in the 
age of the Universe. However, that is not the case for this 
sample. 

The sample selected using the SDSS cuts in the MS are 
not very homogeneous. A source of this inhomogeneity is 
the selection of galaxies in apparent magnitude space. Due 
to the width of the redshift bin within which the appar- 
ent colour cuts are defined, different types of galaxies are 
selected near the edges of the redshift bin including a large 
number of galaxies that have more extended SFHs. Inclusion 
of these galaxies decreases the sensitivity of this method for 
determining H(z). Furthermore, the observing efficiency is 
reduced if only the oldest galaxies are used unless they can 
be pre-selected. 



2.2 Selecting a more homogeneous sample 

Since the age-dating experiment assumes a similar SFH for 
LRGs at different redshifts, we explored other selection cri- 
teria that would yield a more homogeneous set of objects. 
The MS includes rest-frame luminosities, and instead of the 
SDSS selection, we chose to use absolute magnitude cuts 
in the rest-frame to select LRG-like objects. In a real age- 
dating experiment, this would require having a large sample 
of galaxies with multi-band photometry and spectroscopic 
(or high quality photometric) redshifts. For z ~ 0.5, this 
selection can easily be done from the SDSS data (e.g. the 
2SLAQ catalog of Cannon et al. 2006). 

The effect of varying the colour and brightness on the 
normalized distribution of mass- weighted ages of simulated 
galaxies is shown in Fig. [5] for galaxies at z = 0.51. We 
demonstrate that by selecting the brightest galaxies, we are 
able to select galaxies with narrower SFHs; that is, galaxies 
in which more of the stars are formed at earlier times. In 
the lower panel, we show the effect of changing the colour- 
cut, demonstrating that the SFHs are not very sensitive to 
the rest-frame colour-cut we choose. Using these results, we 
decided to use a rest-frame colour-cut of B — V > 0.81 and an 
absolute magnitude cut of My < —23. A brighter magnitude 
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Figure 1. Top: The average star formation histories (star forma- 
tion rate vs. look-back time from today) for LRGs selected from 
the MS using the SDSS cuts in Eisenstein et al. (2001). The solid 
line is for LRGs at z = 0.32, the dotted line is for z = 0.46 and 
the dashed line is for z = 0.51. Bottom: The ratio of the average 
star formation histories in the top panel to the average star for- 
mation history for the z = 0.32 sample. The line-styles represent 
the same redshift bins as for the top panel. 



cut does produce an even more homogeneous sample, but 
results in a dramatic drop in the number of sources (a factor 
of 10 decrease with a 0.5 mag increase). 

In Fig. 13] we plot the average SFHs of galaxies selected 
using our rest-frame cuts at four redshifts (upper panel) as 
well as the ratio of SFR to that obtained at z=0.32 (lower 
panel). It is clear, as compared to Fig. [T] that the galaxies 
selected in this way have more homogeneous SFHs, and the 
small change from low to high redshift is monotonic. To 
demonstrate the variation in SFH for galaxies selected in this 
way at a particular redshift, we plot, in Fig. |4j SFHs of 200 
randomly chosen galaxies from the sample at z = 0.51. In 
the lower panels, we break the 100 galaxies into two samples: 
those that are more than 3<r away from the mean and those 
that are within 3<r of the mean. 

The formation properties of these galaxies are compared 
to the SDSS sample in Table 1. On average, this population 
shows similar dispersion to the objects selected from the 
SDSS cuts. However, a smooth and gradual change is evident 
for our cut as opposed to the dramatic change seen using the 
SDSS cuts. In addition, a constant and lower percentage of 
objects have had significant events of star formation in the 
past as opposed to the SDSS cuts. 

To investigate further the properties of the LRGs se- 
lected with the new cuts, we plot, in the top panel of Fig. [5] 
the mass of the halo which hosts the galaxies vs. the mass- 
weighted age of our selected LRGs for four snapshots. We 
note the more massive halos contain older galaxies, a corre- 
lation which arises naturally in hierarchical structure forma- 
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Table 1. Formation Properties of LRGs 



Cut z Zf z p z e Age mw fsfr^. f3cr 



SDSS 


0.32 


11.86 ± 1.30 


4.95 ± 1.26 


1.46 ± 0.49 


8.32 ±0.28 


0.32 


0.13 


SDSS 


0.46 


11.94 ± 1.23 


5.06 ± 1.46 


1.44 ± 0.44 


7.05 ±0.35 


0.03 


0.19 


SDSS 


0.51 


12.17 ±1.16 


5.05 ± 1.18 


1.49 ±0.47 


6.87 ±0.22 


0.29 


0.17 


Our Cut 


0.32 


11.94 ± 1.23 


4.91 ± 1.14 


1.48 ±0.47 


8.29 ±0.27 


0.14 


0.15 


Our Cut 


0.46 


11.92 ± 1.23 


4.96 ± 1.18 


1.54 ±0.48 


7.20 ±0.29 


0.14 


0.15 


Our Cut 


0.51 


11.89 ± 1.25 


5.00 ± 1.21 


1.58 ±0.47 


6.83 ±0.31 


0.14 


0.15 



a The fraction of galaxies with SFR > 5 M^yr 1 within 1 Gyr of the redshift. 




Figure 2. The effect of varying the brightness cut (top) and 
the color cut (bottom) on the mass-weighted age histogram for 
galaxies selected at z=0.51. 

tion scenarios (de Lucia et al. 2006) and corresponds to the 
observed phenomenon of downsizing (Cowie et al. 1996). 

In addition, these galaxies have large stellar masses with 
values typically above 2 x 10 11 M© with a peak value just 
above 10 12 Mq. There is a strong, lower mass cut-off in 
each of the redshift bins due to the absolute magnitude cut 
in the selection function. The overall stellar mass values are 
in good agreement for stellar masses calculated for LRGs by 
Barber et al. (2007), but their data does not show the strong 
correlation between age and stellar mass for SDSS selected 
LRGs that is seen in the de Lucia models. 

When we substitute metallicity for halo mass, we see 
no correlation. The lower panel of Fig. shows a normal- 
ized distribution of the mass-weighted ages of galaxies from 
our selection. Even at very close redshifts, the ages of the 
galaxies are separated into distinct groups, although a small 
tail of younger galaxies is still present. 

For the rest of the paper, we refer to these galaxies 
as LRGs. Galaxies selected using the old cut will be re- 
ferred to as SDSS LRGs. At four different redshift (z = 
0.32,0.46,0.51,0.56), the total number of LRGs extracted 



Figure 3. Top: The average star formation histories (star forma- 
tion rate vs. look-back time from today) for LRGs selected from 
the MS using our revised absolute magnitude cuts. The solid line 
is for LRGs at z = 0.32, the dotted line is for z = 0.46, the 
dashed line is for z = 0.51 and the dot-dashed line is for z = 0.56. 
Bottom: The ratio of the average star formation histories in the 
top panel to the average star formation history for the z = 0.32 
sample. The line-styles represent the same redshift bins as for the 
top panel. 

from the MS are 1705, 1491, 1448, and 1337 galaxies respec- 
tively. 

2.3 Selecting LRGs from Durham Models 

We applied both the SDSS and the absolute magnitude cuts 
to the Durham model. Although similar density of objects 
were found in the Durham model at z=0.51 by using the 
SDSS cuts, the density of LRGs after the absolute magni- 
tude cut was down by an order of magnitude as compared 
to both the real object density and the de Lucia model. 
This was due to the lack of bright objects in the Durham 
model at intermediate redshifts. This was previously shown 
by Almeida et al. (2008) in the mismatch between the ob- 
served luminosity function of 2SLAQ LRGs and the Durham 
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Figure 4. Top: Star formation histories of 200 galaxies selected 
using the revised cuts from the MS. The black line indicates the 
mean star formation history of 1440 galaxies selected using the 
cuts. The error bars on the black line represent the standard devi- 
ation of the mean. Middle: Galaxies from the sample in the top 
panel whose star formation rates differ by more than 3<r from the 
mean star formation history. Bottom: Galaxies from the sample 
in the top panel whose star formation histories differ by less than 
3cr from the mean star formation history. 
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Figure 5. Top: The masses of the halos hosting LRGs, plotted 
against their mass weighted ages. For all four redshifts, a strong 
correlation is found between age and halo mass. Bottom: His- 
togram of the mass-weighted ages for four redshifts. 



models, where the Durham model has an excess of faint 
LRGs. 

For the selection of objects that were found in the 
Durham model, they had an elevated level of star forma- 
tion, bluer colors, less regular stellar histories, and a large 
spread in the stellar ages. These results are consistent with 
results found by Almeida et al. (2008) and Gonzalez et al. 
(2009) in studies of the Durham model and LRGs. Further 
study is necessary to ascertain the reasons for the differences 
between the models and the observations, and for this rea- 
son, we focus on the results from the de Lucia models for 
the remainder of this paper. 



3 RECOVERING H(Z) FROM MILLENNIUM 
LRGS 

In this section, we explore the uncertainties on H(z) related 
to the extended star formation histories of LRGs and ignore 
difficulties associated with measuring the age of LRGs. We 
use only the galaxies selected with our rest-frame magnitude 
and color cuts. We use the mass-weighted ages, which are 
shown in the bottom panel of Fig. [5] for galaxies in the 
four redshift snapshots initially used. In addition to these 
redshifts, we extract another set of galaxies from snapshots 
with < Z < 1, where we have adjusted our color cuts 
relative to z — 0.51 to account for passive evolution. 

The distribution of the ages, despite refining the cuts to 
provide a narrow distribution, are still highly non-Gaussian 
with a tail reaching to younger ages. With this in mind, we 
compare three different methods to determine the charac- 
teristic age at a given redshift: (1) Calculating the average 
of the ages, (2) fitting a function to the distribution of ages, 
and (3) matching pairs of galaxies. For each approach, we 
calculate H(z) using equation fl}. If we assume no error in 
the redshift, the error in H(z) will only depend on the age 
at the two redshift bins: 



H(zY (i!-t 2 ) 2 



(2) 



In Fig. [6] we present estimates of H (z) according to vari- 
ous methods. We plot the difference between our calculated 
value, tt|, and the expected value for H(z) based on equa- 
tion (3) in Jimenez & Loeb (2002): 



-(i + = -(i + z) 5/2 [n m (o) + n de (o) 

tl 



x exp(3 



dz 



1/2 



(3) 



where Wde = Pde/pde is the equation of state parameter for 
the dark energy, and in a ACDM model of de Lucia et al. 
(2006), w de = -1. 

The average method (top panel) recovers H(z) over a 
range of redshifts reasonable well. If we assume an error in 
the mean age at a specific redshift of 0.03 Gyr (this error 
will be discussed in ^5}, then H(z) can be calculated to a 
precision of 1.6% at z « 0.42 (using the redshift interval be- 
tween z\ — 0.32 and Z2 = 0.51). As can be seen in equation 
J2J, the error in H(z) goes as the inverse of the difference 
in age: the closer the redshift bins, the larger the error. Be- 
tween z = 0.51 and z — 0.56, a small, systematic error in the 
mean age of as little as 0.5% will result in the calculation 
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of H{z) to be off by 10%. Because of the small curvature of 
H(z) over these redshift ranges, large values in Az do not 
introduce a systematic difference in the value of H(z). 

The middle panel presents the results for fitting the 
distribution of ages to a function and determining a char- 
acteristic age. This is akin to fitting the envelope of oldest 
galaxies as in Jimenez et al. (2003). The data were fitted to 
the following function: 



P(t) 



a* g(t) 2 



where, 



(t) = 



t -t 



(4) 



(5) 



For the first test, we allowed a, b, and c to be free parame- 
ters along with t . In the second test, we set b = 5.52 and 
c = 0.66, which were the average parameters found during 
the free fits, and solved for a from the best fit value of t . 
The results for H(z) between the free and fixed fits were 
very similar and Fig. [6] (middle panel) shows the results for 
the fixed fit. Unlike the average estimate, H{z) is slightly 
underestimated at the highest redshifts. This is probably 
a result of the fitting function in equation (4) not fitting 
the normalized age distribution at higher redshifts well. At 
z ~ 0.42, H(z) can be calculated to 1.1% by using our full 
sample at z — 0.32,0.51. 

Finally, we calculate the value of H(z) at a given red- 
shift based on the distribution of H(z) calculated between 
each pair of galaxies in the different redshift bins. If we used 
the entire sample, overlap between the age distributions sig- 
nificantly skews the results. Less than 10 galaxies results 
in anomalous values of H{z). Using the 20 oldest galaxies 
minimizes both the systematic and random errors on the cal- 
culation of H(z), however the results are not substantially 
different for using the 1000 oldest galaxies. Once we have 
calculated the value of H(z) from all of the pairs, H(z) at a 
given redshift is the three-sigma clipped mean of the distri- 
bution. The error is then given by the standard deviation of 
the sigma-clipped distribution. The H(z) estimate at high 
redshifts is also slightly below the model, but well within 
the errors. At z ~ 0.42, H(z) can be calculated to 2.8% 

It is reassuring that we can recover H(z) at a a single 
redshift using various metrics for the characteristic age of 
the population. Even the very gross estimate of the mean of 
the distribution provides very high accuracy (< 2%) on the 
calculation of H(z). Obviously using the ages from the MS 
is a best case scenario where the galaxy ages are determined 
without any error. In the next section, we explore the errors 
associated with age-dating LRGs. 



4 MODELLING AND AGE-DATING LRG 
SPECTRA 

There are a number of stellar population synthesis codes 
available which can be used to generate synthetic spectra 
once a star formation history has been determined (e.g. 
Fioc and Rocca-Volmerange 1997, Bruzual & Chariot 2003) . 
However, the underlying physics is not always well un- 
derstood: thermally-pulsating AGB stars may not be cor- 
rectly modelled which effects infrared magnitudes signif- 
icantly (Maraston 2005); an evolving Initial Mass Func- 



tion (IMF) might be required to model cluster red-sequence 
galaxies which include many LRGs (van Dokkum 2008); re- 
lated to the IMF, the alpha-element enhancement of old 
galaxies is not correctly modelled (Walcher et al. 2009); and, 
finally, dust modelling remains uncertain. Newer generations 
of models are currently being developed and tested (e.g. 
Maraston 2005, Conroy et al. 2009) along with improve- 
ments in the underlying stellar evolution models (Marigo 
et al. 2008). In addition, there is also uncertainty in the 
best method for deriving accurate ages both at low redshift 
(Kannappan & Gawiser 2007, Wolf et al. 2007, Trager & 
Somerville 2009) and at high redshift (Longhetti & Saracco 
2009, Maraston et al. 2009, Muzzin et al. 2009). Unfortu- 
nately, the focus in these studies is more on broad-band and 
low-resolution spectra rather than the high-resolution, high- 
signal to noise (SNR) spectra we expect to need. 

In this paper, we carry out a simple study of the age- 
dating question and leave a detailed discussion for the com- 
panion paper, including a comparison of different population 
synthesis models, age-dating techniques, and comments on 
IMF evolution, alpha-element enhancement and the role of 
AGB stars. 



4.1 Modelling LRG Spectra 

We build synthetic spectra of LRGs by combining single 
stellar population (SSP) spectral libraries of BC2003. The 
libraries are tabulated at ages ranging from 10 s years to 
2 x 10 10 years at a resolution of ~3 A across the whole 
wavelength range from 3200 to 9500 A for a wide range of 
metallicities (Z = 0.0001 to Z — 0.05). Spectral coverage 
over a larger wavelength range, from 91 A to 160 pm is also 
available but at lower resolution. For our studies here we 
have assumed the standard Salpeter initial mass function, 
with mass cut-offs at 0.1 and 100 Mq. We refer the reader 
to BC2003 for details about the stellar evolution prescription 
used in constructing their libraries. 

From the semi-analytic models of galaxy evolution 
based on the MS, one can extract both the star formation 
rate and metallicity of the cold gas out of which the stars are 
formed, as a function of time. Using these and the BC2003 
SSP spectral libraries as input we calculate the emergent 
spectrum according to 



F(X) 



S(t)F SS p(X,t,Z(t))dt. 



(6) 



In this equation, F(X) is the emergent spectrum in the 
rest frame; t, the look-back time from the redshift of in- 
terest; tf orm , the look-back time at the start of formation' 
S(i), the star formation rate per unit mass per unit time; 
Fssp(^,t, Z(t)), the spectrum of an SSP as a function of 
age; and metallicity normalized to unit mass of stars and 
Z(t) is the metallicity as a function of look back time. We 
ignore the effects of dust, since these LRGs are assumed 
to be devoid of any remaining gas and dust (Barber et al. 
2007). 

We generate spectra for all 1705, 1491, 1448, and 1337 
galaxies selected from the MS at z = 0.32,0.46,0.51,0.56. 
In particular, we use the z = 0.51 galaxies to examine the 
limitations in age dating the stellar populations. 

We used the revised colour-magnitude cuts discussed in 
i]2.2l to extract a sample of galaxies from the SDSS catalogue 
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Figure 7. Average LRG spectrum at z=0.46 created by using 
200 galaxies that satisfy the revised cuts in A2.2I from the SDSS 
compared to the best-fit model spectrum produced using BC2003 
with the SFH from the MS. 



(Abazajian et al. 2009). In Fig. [7] we present the average 
spectrum of the first 200 galaxies that match our cuts. For 
comparison, we plot the best fit model spectrum at z = 0.46. 
The fit was made from 3500-6000 A, which includes only the 
high-resolution portion of the models. The best fit model 
has a mass- weighted age of 7.11 Gyrs. The overall shape of 
the model is in good agreement with the observed spectrum 
although there are some inconsistencies at the edges. 



4.2 Age-dating simulated LRG spectra with SSPs 

In our first attempt to age-date simulated LRGs, we deter- 
mine the ages for the 1448 model galaxies at z = 0.51 using 
the SSP template library. For each source, we fit the full 
spectrum from A = 3500 - 9000 A using the SSP library 
described in We do not add any flux errors or noise to 
the spectra and determine the best fit by minimizing the \ 2 ■ 
For our entire sample, the systematic and random er- 
rors from fitting the SSPs are listed in Table 2 as a func- 
tion of redshift. It is not surprising to find a systematic er- 
ror in the age determination using SSPs as compared to 
the mass-weighted age as the most appropriate compari- 
son would be to the light-weighted age. Trager & Somerville 
(2009) perform a similar comparison for red galaxies in the 
Coma cluster comparing observed and simulated popula- 
tions with SSPs. They also find that the SSPs underesti- 
mate the age of the galaxies and the SSP age is poorly 
correlated with the mass-weighted age due to recent star 
formation dominating the spectra. We also find that the off- 
set in the systematic error does not correlate with the age of 
the Universe. The differential ages are less sensitive to sys- 
tematic errors (Jimenez et al. 2004), and if the systematic 
errors were constant, this would not affect the calculation of 
H(z). Since the bias does change with redshift, this prevents 
SSP-determined ages from being very useful for determining 
H(z). 

Based on the sample from McCarthy et al. (2004) of 
early type galaxies, Simon et al. (2005) claim they are well 
fit by a single burst population and use these ages to recover 
a value of H{z). This is consistent with our results presented 
here as the resolution of their spectroscopic data (17 A) and 
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Table 2 


. Errors from the Age Dating 














SSP Fit 






Model Fit 
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Age_^ 


Age 


Ai 




Age 


A b 


(7° 


0.32 


10.08 


8.15 


0.14 


1.85 


8.32 


0.01 


0.29 


0.46 


8.95 


7.42 


0.20 


1.72 


7.22 


0.01 


0.33 


0.51 


8.59 


7.06 


0.23 


1.65 


6.83 


0.02 


0.33 


0.56 


8.25 


6.62 


0.16 


1.69 


6.48 


0.01 


0.33 



a Age of the Universe in Gyrs. 

^ Mean difference between measured age and mass-weighted age. 
c Standard deviation of difference between measured age and 
mass-weighted age. 




SNR 

Figure 8. Top: The average systematic offset between the model 
age and the mass-weighted age as a function of SNR and resolu- 
tion for 10 different galaxies. Bottom: Random error in the age 
determination for the same ten galaxies 



the precision on H(z) are much greater than those that we 
examine here. McCarthy et al. (2004) do fit a majority of 
their objects with burst of less than 0.1 Gyrs, but the objects 
have a range of ages and formation redshifts that vary as 
much as 3 Gyrs. This is consistent with SSP fitting being 
dominated by the most recent star formation. 

To explore the SSP fitting in more detail and to gain in- 
sight into the effects of changing the SNR and resolution of 
spectra, we performed Monte Carlo simulations for a subset 
of 10 galaxies. We selected ten random model galaxy spectra 
from our sample at z — 0.51. The template library and the 
spectra were convolved to a resolution of AA = 3, 5, 10, 20 
A. Noise was added to each spectrum to give a range of sig- 
nal to noise from SNR = 3, 5, 10, 30, 50, 100, 200. The SNR 
is here defined as the average SNR per resolution element 
of the spectrum between 3000 and 9000 A, assuming simple 
shot noise. \ 2 minimization was used to find the metallic- 
ity, age and normalization of the best fitting SSP spectrum. 
For each galaxy, this process was repeated 1000 times. Fig. 
[8] shows the average of the results from these simulations. 
As an example of the typical spread in ages, at a value of 
SNR — 100 and AA = 3 A, the minimum and maximum 
values for the ten galaxies were A age = 0.148 — 1.23 Gyrs 
and a age = 0.025 - 0.045 Gyrs . 

The precision in this method, even when using the SSP, 
improves with increasing SNR as expected, and this method 
can very reliably reproduce the same age at high SNR for 
all resolutions. The random error on the age does decrease 
for smaller resolutions allowing for a better determination 
of the age. However, the poor accuracy of SSPs is again 
highlighted in Fig. [8] The mean offset in age seems to settle 
at a given value at SNR > 30, but it does not seem to be a 
monotonic function of the spectral resolution. 

The accuracy of the age-dating is clearly critical for cal- 
culating H(z). Unless systematic biases are constant with 
redshift and can be subtracted out or are well behaved and 
can be removed with modelling, they will contribute a sig- 
nificant systematic error to the age calculation. Even in this 
relatively straight forward simulation of using model galax- 
ies from the MS, the SSPs are not able to reproduce the age 
without significant bias and are thus not adequate for the 
cosmic chronometer method. 



4.3 Age-dating LRGs with model spectra 

The large model library of different LRGs that we have 
created provide a mapping from spectra to mass-weighted 
ages. We can potentially use these spectra as our templates 
instead of the SSP templates. To test this, we extracted 
one spectrum from the models and then used the remain- 
ing models to age date that spectrum. Then, we repeated 
this for all of the other spectra. For each spectrum, we used 
its mass weighted age as its fiducial age. Testing the entire 
sample for 1448 z = 0.51 model galaxies, we find a differ- 
ence between the mass-weighted age and the measured age 
of A mm = —0.02 Gyr with a dispersion of a mm = 0.32 Gyrs. 
The result for the other redshifts are listed in Table 2. The 
systematic error, A mw , and dispersion, a mw , have been sub- 
stantially reduced as compared to using SSPs. In addition, 
these ages can be used to measure H(z) because they have 
very small, regular bias with respect to the age of the Uni- 
verse. 



5 AN OBSERVING PROGRAM 

In this section, we explore the minimum observing time re- 
quired to recover H(z) to a precision of 3,5, and 10% at 
2 ~ 0.42. Table 2 indicates that uncertainties on individ- 
ual ages of galaxies could be as low as 0.3 Gyr if suit- 
able templates can be used. However, since it is still un- 
clear what the uncertainty on individual ages of galaxies 
will be in a realistic experiment, we explore observing re- 
quirements for four values of this uncertainty (0.05 Gyr, 0.5 
Gyr, 1 Gyr and 2 Gyr). We consider observations at two 
redshifts: z — 0.32 and z = 0.51, giving a redshift interval of 
Az — 0.19. To simplify our estimate, the mean ages at the 
two redshifts are calculated using the average age of LRGs 
at each redshift (the first method discussed in §3). Accord- 
ing to equation ([2]), the uncertainty on the mean ages will 
have to be <J <age> — 0.03,0.05,0.10 Gyr to measure H(z) 
with 3, 5, 10% precision. 

To estimate the uncertainty in mean age as a function 
of the number of galaxies and the uncertainty on individual 
ages, we do a simple Monte Carlo simulation. We assume 
the uncertainty on individual galaxy ages are normally dis- 
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Figure 9. The uncertainty in the mean age (in Gyr) of LRGs 
at z = 0.51 as a function of the number of galaxies used for the 
measurement. Wc plot curves assuming four different uncertain- 
ties (see legend) in the measured age of an individual galaxy in 
the sample. 



Figure 10. The total observing time required with SALT to mea- 
sure H(z) at z = 0.415 to 3, 5, and 10% as a function of signal to 
noise of the observations. All observations are for AA = 3A and 
an overhead of 300 s per observation. 



tributed, and the galaxy ages are drawn from the probabil- 
ity distribution for galaxies at z — 0.51 given by the fit to 
equation Q to the normalized age distribution. For each 
N, the simulation is repeated 1000 times and the standard 
deviation in the mean age is calculated. The results are pre- 
sented in Fig. and are used to calculate the total number 
of galaxies required to reach our desired precision. 

The other constraint on the total observing time is the 
exposure time per galaxy. We assume that we are only able 
to measure one galaxy at a time but note that multi-object 
spectroscopy could be used. In Fig. [5] the SNR for different 
spectral resolutions was presented. Even though the results 
in Fig.[8]are obtained using SSP fitting, we use them to pro- 
vide a crude estimate of the resolution and signal to noise 
required to derive the random error on the ages of individ- 
ual galaxies. In Fig. 1101 we plot the total time needed as a 
function of signal to noise to calculate H(z) with different 
precision using the Robert Stobie Spectrograph (RSS, Kob- 
ulnicky et al. 2003) on the Southern African Large Telescope. 
We estimate that a 10% measurement is feasible using less 
than 20 hours of observing time but that a 3% measurement 
could require ~ 180 hours, even with our fairly optimistic 
estimates of a ag s- 

We note however, that the RSS can be used in multi- 
object mode and, depending on the clustering, several 
LRGs may be expected in each field of view. Using our 
cuts, we estimate the space density of LRGs to be 3.9 x 
10" 5 Mpc~ 3 from the SDSS (without considering the ef- 
fects of incompleteness), which is comparable to the value 
of 3.5 x 10~ 5 Mpc~ 3 from the MS. If we are able to observe 
two additional LRGs per set-up with z — 0.1 — 0.6 while 
calculating H(z) to 3% as outlined above, it would give us 
a sufficient numbers of LRGs to calculate H(z) to 10% at 
redshifts between z = 0.1 — 0.6. This would put a far tighter 
constraint on the value of H(z) than measuring it from two 
redshift bins alone. 

6 CONCLUSIONS 

In this paper, we have shown that within the self-consistent 
Universe of the Millennium Simulation that Luminous Red 



Galaxies, when selected in an appropriate manner, can be 
used as cosmic chronometers. By selecting the galaxies from 
their rest-frame properties, we find that we can create a 
more homogeneous sample of objects than the apparent 
magnitude cuts used by Eisenstein et al. (2001). For galax- 
ies at z ~ 0.3 — 0.6, we find that cuts of My < —23 and 
(B — V)o > 0.81 select a sample of galaxies with similar star 
formation histories and formation redshifts in the MS. 

The galaxies selected in our final cut do show very sim- 
ilar star formation histories with very few of the galaxies 
showing any star formation since z ~ 1.7. The ages also 
show something akin to downsizing (Cowie et al. 1996): the 
oldest galaxies are in the most massive halos. The distribu- 
tion of the ages are very similar with each showing a small 
tail towards younger age, but with the very strong peak at 
a single age. 

However we do note limitations in the use of the models 
of de Lucia et al. (2006). They are able to reproduce some 
of the major trends seen in local galaxy populations and 
seem to have similar number density as the SDSS out to 
z ~ 0.5. However, the model does have some limitations in 
reproducing some of the observed properties of LRGs while 
alternative models produce equally acceptable fits to the 
data. One such model explored here (Bower et al. 2006) was 
unsatisfactory in matching the observed number density of 
LRGs defined by an absolute magnitude cut at a z ~ 0.5. 

Even though these galaxies have extended SFHs, they 
can be used as cosmic chronometers to recover the cosmol- 
ogy used in the MS. Using only galaxies selected from two 
redshifts, H{z) can be calculated to a precision of less than 
3% using three different methods. All three methods (av- 
eraging the ages, calculating a fit to the distribution, and 
comparing the ages of pairs of galaxies) were able to recover 
the cosmology used in de Lucia et al. (2006). 

In Sj4] we showed that SSPs were not sufficient to ac- 
curately recover the ages for individual galaxies, which has 
also been shown for other samples (Maraston et al. 2009, 
Trager & Somerville 2009) . Despite the relatively simple na- 
ture of the SFHs in Fig. |3j the calculated SSP ages were 
dominated by the most recent burst of star formation. If we 
use the average star formation history from the MS, we are 
able to replicate the properties of the model spectra, but in 
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the next paper, we look in far greater detail at the modelling 
and fitting of LRG spectra. However, the semi-analytic mod- 
els likely overestimate the extent of star formation in LRGs, 
and SSP models may provide better fits than indicated here. 

Finally, we estimated the required time to complete the 
project using RSS on SALT. If systematics in the age can be 
controlled, H(z) at z ~ 0.42 can be calculated to 3% in a to- 
tal of ~ 180 hours, which includes observing overheads. It is 
likely that tighter constraints could be put on the evolution 
of H(z) with additional data that could be obtained while 
making these observations. In addition to constraining H(z), 
it would contribute a wealth of detailed information about 
the evolution of the most massive galaxies at intermediate 
redshift. 

Throughout this work, we have highlighted a number 
of assumptions that we have made. Even with our refined 
selection, the star formation histories of these galaxies are 
not perfectly homogeneous as assumed by Jimenez & Loeb 
(2002), but H(z) can still be calculated to an accuracy better 
than 3%. Other improvements on the measurement can be 
made by selecting a larger or more homogeneous sample 
as done by Stern et al. (2010). Uncertainties in the star 
formation history of simulated galaxies introduce far smaller 
errors than the contribution associated with age-dating the 
galaxies. In our next paper, we explore age-dating in more 
detail and note that increased uncertainties resulting from 
age-dating may lead to increased observing requirements. 
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